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Relativistic Weyl fermion (WF) often appears in the band structure of three dimensional magnetic 
materials and acts as a source or sink of the Berry curvature, i.e., the (anti-) monop ole. It has been 
believed that the WFs are stable due to their topological indices except when two Weyl fermions 
of opposite chiralities annihilate pairwise. Here, we theoretically show for a model including the 
electron-electron interaction that the Mott gap opens for each WF without violating the topological 
stability, leading to a topological Mott insulator dubbed Weyl Mott insulator (WMI). This WMI is 
characterized by several novel features such as (i) energy gaps in the angle-resolved photo-emission 
spectroscopy (ARPES) and the optical conductivity, (ii) the nonvanishing Hall conductance, and (hi) 
the Fermi arc on the surface with the penetration depth diverging as approaching to the momentum 
at which the Weyl point is projected. Experimental detection of the WMI by distinguishing from 
conventional Mott insulators is discussed with possible relevance to pyrochlore iridates. 


Introduction — Weyl fermions (WFs) in solids attract 
recent intensive interests from the viewpoint of their 
novel quantum transport properties and chiral anomaly. 
The WF is described by the 2-component spinors orig¬ 
inating from 4-component Dirac spinor when the mass 
m is zero in the Dirac equation. The realization of WFs 
in condensed matters has been recently established [1- 
3]. In magnetic materials, the time-reversal symmetry 
is broken and the energy dispersion of Bloch wavefunc- 
tion has no Kramer’s degeneracy. In this case, the band 
crossings between the two bands are described by a 2 x 
2 Hamiltonian as 

H{k)=eo{k)T X (1) 

i=l,2,3 

with 2x2 Pauli matrices a* {i = 1,2,3). Three condi¬ 
tions of the band crossing hi{k) = 0 for {i = 1,2,3) can 
be satisfied in general by appropriately choosing the three 
components of the crystal momentum k. Weyl points 
sometimes exist exactly at the Fermi energy when dic¬ 
tated by some symmetry and topology of the Bloch wave- 
functions, for example, in Dirac semimetals [4-8]. More 
recently, experimental discovery of Weyl semimetals in an 
inversion broken material TaAs has been reported [9-11]. 

Weyl fermion plays an important role in the con¬ 
text of the Berry phase, which is defined by = 

—'i{unk\^k\unk) (I'^nfc) • the periodic part of the Bloch 
wave function with the band index n = ± and the mo¬ 
mentum k ) and acts as the vector potential in the mo¬ 
mentum space. The Berry curvature b^k = x a^k is 
the emergent magnetic field, and can be enhanced near 
the band crossing points. When one expand Eq. (I) 
around the band crossing point (Weyl point, which we 
assume to be k^ = 0), there appears the WF described 
with hi(k) = rjUpki, by an appropriate choice of the co¬ 
ordinate /ci’s, where vp is the Fermi velocity. The sign 
T] = ±1 specifies the chirality of the WF, and the Berry 
curvature of the lower eigenstates (n = —) of the Hamil¬ 


tonian in Eq. (I) is obtained as 


which diverges as \k\ 0 and the total flux pene¬ 

trating the surface S enclosing the Weyl point is given 
by dS ' b-k = This indicates that the WF 

acts as the magnetic monopole (anti-monopole) for = I 
(77 = —I); the magnetic charge ^ plays a role of 

topological index. Strong Berry curvature leads to the 
enhanced anomalous Hall effect [I] as well as the chiral 
magnetic effect which results in the negative magneto¬ 
resistance [12]. 

Figure I shows the schematic figure of the three dimen¬ 
sional first Brillouin zone in which two WFs exist along 
the kz direction. One can define the Chern number 

CKK) = I (3) 

for the plane of fixed kz. When we consider Ch{kz) as a 
function of kz, there appears the jump by rj dXkz = 
i.e., -coordinate of the Weyl points. Therefore, due to 
the periodicity of Ch{kz) by kz kz P ‘1 'k jc, we need 
the pair of 77 = I and —I [13, 14]. The existence of a 
single (an odd number of) WF is also excluded. There¬ 
fore, the annihilation of a single WF is prohibited, i.e., 
the only way to destroy the WFs is to annihilate a pair of 
WFs with opposite chiralities either by making the two 
WFs approach to each other in the momentum space or 
by introducing a scattering between two WFs with some 
density-wave-type order. The former scenario is actually 
proposed for the transition between the Weyl metal and 
insulator in pyrochlore compounds [15]. The latter one 
is also discussed intensively [16, 17]. Meanwhile, effects 
of the electron correlation have been discussed for WFs 
by several methods including random phase approxima¬ 
tion [18] and, more recently, cluster perturbation the¬ 
ory [19]. However, the possibility of the Mott gap open¬ 
ing at each WF has never been explored thus far to the 
best of the present authors’ knowledge. 
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FIG. 1. Schematic picture for the Weyl semimetal. A Weyl 
point plays a role of a source or sink of the Berry curvature, 
i.e., the (anti-)monopole in the momentum space. A pair 
of Weyl points with opposite charges is accompanied with a 
Fermi arc (the curve on the right side surface). 


In this paper, we study the effect of the electron cor¬ 
relation U on WFs by using a simple model which is 
exactly solvable. It is shown that the Mott gap due to 
U open at each WF without the pair annihilation, while 
the topological properties are kept unchanged. Namely, 
the magnetic charge of the WF is unchanged with the 
role of poles in Green’s function being replaced by its ze¬ 
ros. The Hall conductance (Jxy remains nonvanishing and 
the Fermi arc on the surface remains, while the Green’s 
function and the optical conductivity (7xx{^) show the 
gap. Therefore, this Mott insulating state is identified 
as a topological Mott insulator, and we name it “Weyl 
Mott Insulator (WMI)”. The experimental detection of 
this new state is also discussed. 

Model and Green’s funetion — The model we study is 
given by the Hamiltonian 


^ = E 


ipiHk) ■ a'lp + Ic(nfc - 1)^ 


( 4 ) 


where is the two-component spinor, 

and Uk = The most peculiar 

nature of this model arises from the electron-electron 
interaction which is local in k, i.e., the Hamiltonian 
is decomposed into independent /c-sectors. In the real 
space, this corresponds to the non-local interaction in 
the limit of forward scattering. A similar idea has been 
explored to study the Mott transition [20] and the spin- 
charge separation [21]. This locality of the interaction 
in k enables the exact solution of this problem. One 
can introduce the unitary transformation U{k) satisfy¬ 
ing U{ky[h{k) • a]U{k) = cr^h{k) with h{k) = \h{k)\ 


(a) \ r // \ / (b) 
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FIG. 2. Energy spectrum of WMI. (a) Poles of the Green’s 
function. Dashed line represents the band structure of a non¬ 
interacting Weyl semimetal. Energy bands are shifted by 
±f//2 and show a Mott gap of U. (b) Excitation spectrum 
indicating the region of non-zero conductivity a(q, cu). Excita¬ 
tion gap of H + is finite for nonzero momentum transfer 
q. The spectrum is singular at g = 0 where gapless excitations 
are allowed due to /c-local excitations without the energy cost 
of U. 


and a new spinor (j)^ = [/(/c)Vfc = and 

then obtain 


hk{nk+-nk-) + ]^U{nu-lf , 

( 5 ) 

with nk± = a\^^ak± and = nk+ + Uk-- There are 
four eigenstates and eigenenergies: (i) |vac) with E = 
(ii) aj[^|vac) with E = h{k)^ (hi) aj[._|vac) with 

E = —h{k), and (iv) |vac) with E = 

Using these solutions, one can easily obtain the ther¬ 
mal Green’s function in the zero temperature limit as 

G~^{k, iuj) = iujl — /left(^) • cr (6) 


U{k)^HU{k) = ^ 


where hQf^{k) = n{k)[h{k) + 17/2] with n{k) = 

h{k)/h{k). (For details, see Supplementary Information 
SI.) As can be seen from Eq. (6), the energy dispersions 
of the poles are e±{k) = Eh{k) + U/2, where the Mott 
gap of U exists even at the Weyl point with h{k) = 0 as 
shown in Fig. 2(a), which can be measured in the angle- 
resolved photoemission spectroscopy (ARPES). It should 
be noticed that Eq. (6) is derived from the exact Green’s 
function and is not a result of some mean-field approx¬ 
imation. Thus, the WFs disappear due to the electron 
correlation without the pair annihilation. 

Topologieal properties — The topological index for the 
interacting electronic systems can be defined in terms of 
Green’s function [22]. For the (2+l)D case, it is given by 



X t^[{G-^dkA{G-^dk,G){G-^dkM^ ( 7 ) 

where run over 0,1, 2, and £ai3'y is the totally an¬ 

tisymmetric tensor. Plugging Eq. (6) into Eq. (7), one 
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FIG. 3. Temperature dependence of the optical conductiv¬ 
ity of WMIs. We plotted results for ^ = 10 (blue), /3 = 3 
(green), and (3 — 2 (red). The delta functions at cj = 0 and U 
from the intraband contributions are denoted by bars whose 
heights represent mutual ratios of the weights. The dashed 
line represents cr{uj) for a free WF. We set = 1. 


obtains 


Ch{k,) = 


f 

J 


d^k 



dnk\ 

4j’ 


(8) 


where the /c-integral is over kx and ky for fixed kz. Since 
n{k) = hefi{k)/\hefi{k)\ = h{k)/\h{k)\, Ch{kz) does not 
change in spite of the gap opening at Weyl points. The 
Green’s function in Eq. (6) depends on the direction in 
which k approaches to the Weyl point ko and still plays a 
role of a source (sink) of the Berry curvature. At exactly 
/co, G has a zero at cj = 0 when one averages over the 
direction of limitk^h^n{k). Namely, the role of a pole is 
replaced by a zero in the topological properties of Green’s 
function [22]. Because of the bulk-edge correspondence, 
nonzero topological index Ch{kz) indicates that the ex¬ 
istence of the surface states on the side surface, i.e., the 
Fermi arc. (For details, see Supplementary Information 
SIT) 

Therefore, the present insulating state is topological 
and we call it “Weyl Mott insulator (WMI)” distinct from 
the usual antiferromagnetic Mott insulator (AFI). 

Optical conductivity — Now we study the conductiv¬ 
ity, which is given by the two-particle correlation func¬ 
tion. We consider a single WF with the Fermi velocity 
vf described by h{k) = vpk. It is crucial to distin¬ 
guish between the nonzero momentum transfer q and 
exactly q = 0. In the former case, the double occu¬ 
pancy of the electrons will be created, while not in the 
latter case, which brings about the singularity or discon¬ 
tinuity at gr = 0. This reflects the long-range nature 
of our Coulomb interaction in Eq. (4). Eor gr ^ 0 , the 
particle-hole continuum starts from uj = U VF\q\ for 
the transition of an electron from k to k q. 


(a) (b) 
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FIG. 4. Phase diagrams of the Weyl semimetal, (a) The phase 
diagram of the Hamiltonian in Eq. (4) for WMIs. (b) The 
conjectured phase diagram for Weyl semimetals with more 
realistic interactions. Insets show the excitation spectra of 
the conductivity a(q,(ju). The singularity at g = 0 in the 
excitation spectrum for the WMI phase turns into a gapless 
collective excitation for realistic interactions. 


Eor the optical conductivity, the momentum of the in¬ 
cident light q is finite although small, and hence we take 
the limit gr ^ 0. In this limit, the optical conductivity 
at zero temperature for a single WE is obtained as (See 
Supplementary Information Sill for details) 

where 0{x) = l(x > 0) and 0{x) = 0{x < 0). The optical 
conductivity shows a Mott gap of U and its asymptotic 
behavior for large uj is given by a{uj) ^ jVlhvF^aj 
which coincides with the well known result for a free WE 
[23]. In Eig. 3, the optical conductivity is plotted for 
various temperatures (See Supplementary Information 
Sill for the derivation). As the temperature increases, 
peaks at cj = 0 and uj = U appear, as denoted by bars 
whose heights represent mutual ratios of the weights of 
the peaks. The peak at cj = 0 is a Drude peak for finite 
temperatures, while the peak dit uo = U arises from an 
intraband contribution in which a WE at k scattered to 
k-\-q within the same band feels a Coulomb repulsion U. 
The appearance of the peak at uo = U in a{uj) for WMIs 
contrasts to the absence of such a peak for AEIs, because 
the peak at uj = U originates from the correlation effect. 
In addition, the appearance of the in-gap absorption in¬ 
dicates the fragile nature of the Mott gap compared with 
the single-particle band gap. 

Eor the case of exactly g^ = 0, only the vertical transi¬ 
tions within the same /c-sector contribute to the conduc¬ 
tivity as indicated by the red line in Eig. 2(b). Since no 
double occupancy is created in this case, no energy cost 
of U occurs. Therefore, the conductivity cr(gr = 0,ci;) is 
given by Eq. (9) with [7 = 0. 

Discussion — Now the relevance of the present re¬ 
sults to realistic systems is discussed. There are clear 
differences between the WMI and the AEI due to their 
topological nature: (i) The Hall conductance is finite in 
WMI while it is zero in AEI. (ii) Correspondingly, the 
Eermi arc on the surface remains in WMI while not in 
AEI. Also importantly, the phase transition between the 
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WMI and AFI is possible once the former exists as ex¬ 
plained below. Figure 4(a) shows the phase diagram of 
the present model. The horizontal axis is the separation 
Akz between the two WFs which is controlled by, e.g., 
the strength of the antiferromagnetic long range order 
parameter M. When [/ = 0, the phase transition occurs 
from the Weyl semimetal to the AFI by the pair anni¬ 
hilation of WFs at M = Me. Once the interaction U is 
switched on, we always opens the gap and the system be¬ 
comes the WMI as long as Akz is finite. Along the phase 
transition line = 0 at M = Mc^ 1/ > 0, the pair 
annihilation of the two zeros of the Green’s function oc¬ 
curs, which is distinct from that at 1/ = 0 where the two 
poles collide and pair-annihilate. Here one must consider 
the peculiarity of the present model. The effect of the 
long range Coulomb interactions is marginally irrelevant 
as in the case of quantum electrodynamics (QED) [24]. 
As for the short-range Coulomb interaction I/, it is ir¬ 
relevant. This means that there must be a finite range 
of Coulomb interaction within which the WFs remains 
gapless and stable. On the other hand, the strong U 
limit in the lattice model corresponds to the localized 
electron and hence a trivial AFI. Therefore, the conjec¬ 
tured phase diagram of a more realistic model is given in 
Fig. 4(b), where the successive transitions from the Weyl 
semimetal to the WMI, and from the WMI to the AFI 
occurs as the strength of the interaction increases. (The 
separation of two WFs is reduced also as the interaction 
increases and hence the trajectory should goes as U and 
M simultaneously increase.) 

Now we discuss the Green’s function and the two- 
particle correlation function for realistic electron-electron 
interactions given by 

Hc= XI (10) 

The self-energy T^{k = 0,co’) of the Green’s function in 
the second order in V{q) is given in Supplementary In¬ 
formation SIV. It is concluded that the gap of the spectral 
function is stable and remains nonzero. The two-particle 
correlation functions such as a{q,uj), on the other hand, 
is gapless at gr = 0 for the Hamiltonian in Eq. (4). Eor 
a finite size system of N sites, the number of poles form¬ 
ing this gapless excitation in the two particle correlation 
function [the red line in Eig. 2(b)] is of 0{N) (which is 


the number of k where the excitation can be created). 
In general, the number of poles for a collective excitation 
is of 0{N), while that of a continuous excitation arising 
as a pair of single-particle excitations is of 0{N^). This 
infers that the gapless excitations at q = 0 corresponds 
to a collective excitation. Eor realistic interactions, the 
discontinuity at gf = 0 should be removed. In this case, 
we conjecture that the vertical transition [red line in the 
inset of Eig. 4(a)] turns into a collective mode with a lin¬ 
ear dispersion as shown in the blue line in the inset in 
the WMI phase of Eig. 4(b). 

These considerations offer a different scenario to inter¬ 
pret the phase diagram of pyrochlore iridates R 2 IT 2 O 7 [3, 
15]. As the radius of the rare-earth ion R is reduced, the 
correlation strength increases. A recent optical measure¬ 
ment in Nd 2 lr 207 has revealed the opening of the Mott 
gap of the order of 0.05eV [25]. A transport experiment 
also discovered the metallic domain wall states even in 
the Mott insulating state, i.e., the bulk is insulating while 
the domain wall is metallic [26]. As the correlation is fur¬ 
ther reduced, these surface metallic states also disappear. 
One scenario is proposed by Yamaji et al. [27] based on 
the mean field theory. As an alternative scenario, one can 
consider the two types of Mott insulators, i.e., the WMI 
and the AEI, and the disappearance of the metallic do¬ 
main wall states signals the phase transition between the 
two phases. Namely, since two domains of the antiferro¬ 
magnetic order correspond to opposite signs of axy and 
hence the two-dimensional chiral surface modes are ex¬ 
pected to appear at the domain boundary in the WMI 
phase. However, we note that this requires a symmetry 
lowering to violate the cancellation of the Chern vectors 
pointing toward four momentum directions equivalent to 
(1,1,1) which makes dxy zero in the cubic symmetric case 
[28]. The smoking-gun experiment for the WMI should 
be the ARPES to detect the Eermi arc even in the Mott 
insulating phase as mentioned above. 
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Supplementary Information for ^^Weyl Mott insulator” 


SI. GREEN’S FUNCTION 


We derive the Green’s function for the Hamiltonian given by 

H = 'ipl,h{k) ■ (Txpj^ + ^U{nf. - - 1), 


(SI) 


where 



(S2) 


a are Pauli matrices acting on spin degrees of freedom, nj^ is the density operator nj^ = U is the magnitude 

of the repulsive interaction. The repulsive interaction is infinite-ranged in the real space, which can be represented 
in a local way in the momentum space. Thanks to the locality in the momentum space, the Green’s function can be 
exactly computed for this Hamiltonian as follows. 

First we perform a unitary transformation 






that diagonalizes the single-particle part of the Hamiltonian as 

U\k)[h{k) ■ cr]U{k) = h{k)a^. 
Then the Green’s function is transformed as 

G„Mfe,r) = -(T,cfeyr)4p = 

where denotes the time ordering and . 

In this basis, the Hamiltonian is diagonalized as 


where |0) is the vacuum state and h{k) = \h{k)\. Thus the expectation value is written as 

+2e-^^. 


(S3) 


(S4) 


(S5) 


II 

to| q 

(S6) 

HbljO) = -hik)bl_\0), 

(S7) 

Hbl^\0) = h{k)bl^\0), 

(S8) 


(S9) 


(510) 

(511) 


In the right hand side of the equation for ^), the forth term vanishes and other terms are nonzero when 
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a = a'. Thus we obtain 

_ jL(^gr(-/i+f )-/3f )+/3/i^^ 

z 

(&fc_(T)&fe_) = ^((0|&fe_('r)6)^_|0)e“^¥ + (0|6fc_^6/._(T)6^_6^_^|0)e“^'^(^)) 

= JL(^g^(^+¥)-/5¥ 

z 

Therefore, the imaginary-time Green’s function is given by 

G++ik,iuJn) = 

_ 1 + e^'^ + e-^¥ \ 

^ y icOn - h + ^ iw„ - ^ § y ’ 

G__(fc,za;„) = -^^dre*“""(6fc_(r)6j^_) 

_ 1 / e-^'^ + e-^7 e^'^ + e-^¥ \ 

^ I + h — ^ iiUfi ~\~ h -\- ^ j 


A. Green’s function for T = 0 


In the zero temperature {/3 oo), the above Green’ function reduces to 

~ T 


G - (kG^n) 

In the original basis, the Green’s function is given by 
G{k^ iu)fi^ — 7 


i. u ’ 

lUJrn. — a — ^ 


iujfi h ^ 


iujn — (h + n • (T ’ 


h 

^=h- 


(S12) 


(S13) 


(S14) 


(S15) 


(516) 

(517) 


(S18) 


SII. FERMI ARC 


In this section, we study the Fermi arc in the WMIs. Nonvanishing topological indices for the WMIs indicate that 
the Fermi arc remains in the WMIs, which we verify in the following. Since our model [Eq. (SI)] is diagonalized at 
each /c-point and hence we can consider an effective two-dimensional model for each -sector when the system is 
periodic in the ^-direction. Because of the bulk-boundary correspondence, we expect that “edge channels” for each 
kz form a Fermi arc. More explicitly, one can obtain the surface bound state from the effective Hamiltonian 




h{k) + ^ 


n{k) • (T, 


h{k) = VF{kx, ky, kz), (S19) 


by replacing the momenta kx,ky with the derivatives —idx,—idy. Away from the plane kz = ^koz, the surface 
state is almost unchanged from the noninteracting case. The nontrivial issue is how the surface state behaves as kz 
approaches Specifically, the problem is whether the penetration depth of the surface states diverges or not with 

kz ^koz- Intuitively, the finite gap U indicates that the length scale ^ remains finite, i.e., ^ = Hvf/U. However, 
it turns out not when one studies the effective Hamiltonian in Eq. (S19) and the asymptotic behavior of the surface 
bound state as \x\ oo (here we assume ky = 0) by tentatively taking the limit of \kx\ \kz\. In this limit, 
= {vf + ^\kz\~^)[—idxcr^ + kzcr^], which indicates that the penetration depth diverges with ^ = \kz\~^. In any 
case, the length scale is determined by \kz\~^ even when we take into account of the higher orders in dx- Therefore, 
the surface bound states penetrate into the bulk as kz approaches to 



sill. OPTICAL CONDUCTIVITY 


We study the optical conductivity a(cu) for a single WF described h(k) = vpk. In the following, we set the Fermi 
velocity Ui? = 1, which can be always restored by the dimension analysis. 


A. Matrix elements 

Here we calculate matrix elements that we will need in evaluation of conductivities, i.e., (±|cri|±). We first param¬ 
eterize the direction of the momentum as 


n = (sin 0 cos 0, sin 9 cos 0, cos 0). 

Then the wave functions that diagonalize the Hamiltonian are written as 



The matrix elements are given by 


(S20) 


(S21) 


(+ +) = sin ^ COS 0, 

(S22) 

(““) = — sin^cos^. 

(S23) 

= cos6>cos0 -h isin0. 

(S24) 

(+ cTy +) = sin6>sin0. 

(S25) 

(+ cTy —) = —i cos 0 + cos ^ sin 0. 

(S26) 

In the evaluation of the optical conductivity, we need 

f 47r 

J sin9d9d4){±\ax\±){±\cr^\±) = —, 

(S27) 

J sm9d9d(j){+\ax\-){-\(Jx\+) = 

(S28) 

In the evaluation of the Hall conductivity as a function of we need 

J d(l){-\-\ax\-\-){-\-\cry\-\-) — 0 , 

(S29) 

r ^ 

/ d(l){-\-\ax\—){—\cry\-\-) = 27TicosO = 27Ti-^. 

(S30) 


B. Zero temperature 

We first focus on the conductivity a{uj) for the zero temperature. The Green’s function is given by 


G {iujjYi) — ^ 


1 


{ioo^Y-{k + ^Y 
with n = k/\k\. The optical conductivity is given by 

Q{u ie) 


tUJm + I /c + — 1 n • O- 


r(co’) = Re 


— lUJ 

d^k 


/ d'^k 

—— ^ ^ \)T^G ^ jYi^ (7 xG {k -j- “1“ 


(531) 

(532) 

(533) 
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The integrand of Q{iVt) reads 

^ ^ i(jjfYi^(JxG{k^ iidm T 

iuJrn 

= E 


tr 


{itOm + - {k+ ^)2 (iWm)2 - (fc + ^)2 


U 


itOm +i^ + \ k A - ] n ■ a ] I iujm + \ k A - ] n ■ a ] a. 


U 


iuJm 

= E 


{iuJrn + - {k-\- f )2 (icJ^)2 - (/c + 


{iuJm + in)iuJm + ( 1 + ^ ) 


{iuJm + - {k+ f )2 {iuJrn)‘^ - (/c + f )2 


{iu)m T i^^iidrn o ( ^ o 


By using the formula 


^ I" - 7 /, 


[(icj^ + iVt)iujm — (^kc] 


(a + 6)(1 - c) 


[{iuJrn + ~ 2 [if2 — (a + 6)] [if2 + (a + 6)] 

for a > 0, 6 > 0, we perform the summation over iujm foi* the above equation and obtain 

5 ;;tr[ 0 (fc,»„KG(fc,i„„ + i!iK| = ^ ,isi-( 2 t+E» + ( 2 t + G)] ' 


(S34) 


(S35) 


(S36) 


After the analytic continuation iff —> w + ie, only the pole at fc = w — ^ contributes to the imaginary part of the 
A:-integral. Thus, we obtain 


4 C°° ^ 

Im[(3(w)] = J k^dk 


(fc+f) 

LU T (2ik T t/) 


Im 


uj ie — {2k -\-U) \ 


1 




(S37) 


where we used the formula Im 


k—a—ie 


^ ^ = nS(a). Hence, the optical conductivity for the zero temperature is given by 

Im[Q(w)] 


a{Lo) = -- 


OJ 


^ (w - uf 0 {lo-U). 


24:7ruj 


By restoring the unit of jh and the Fermi velocity end up with 

= 


(S38) 


(S39) 


1. Poles ofa{q,(jj) 

Let us study the locus of the poles of the two-particle correlation function that contribute to the conductivity cr(gr, uj) 
for nonzero q. From Eq. (S35) and setting a = \k\^ and b = \k q\ the poles of ^A^{^^'^^rn)o'xG{k + 

gr, iuJrn can be read off as cj = a + 6 = \k\ P \kq\U. By using the formula \k\ P\k-\-q\ > |gr| and restoring 

the Fermi velocity vp, the lower bound of the poles is given by 

uj = U P vfIqI- (S40) 


C. Finite temperature 

In this section, we calculate the optical conductivity a{uj) in the finite temperature. In doing so, we consider 
contributions from interband and intraband transitions separately as 


a{uj) = cr^^^^^{uj) P 


(S41) 
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1. Interhand transition 


The interband contribution to is given by 


Winter 


1 

(27r)^ 


/ 


k'^dkAinterii^), 


(S42) 


where 


Winter 

= J sin 9d9d(l)y^}G++{k,iLOm){+\o'x\-)G — {k,iujm + i^){-\crx\+) 


= ^z-^ 

3 




+ + '^^){-\-\o'x\—)G - {k^ '^^m) {—\o'x\-^)] 

f g^/c 0-/5^ g-/3/c \ _ou_ ( 


s=±l 


‘Zk U isl^ 


+ 


2k + isQ 


e ^ 2 


+ 


e + e ^2 


.-/3/c 


2k + isQ 2k — U isQ 




e + e ^2 




g/3/c g /3 ^ 


e + e ^2 


p-/3¥ 


— (2k T t/) T — 2k T isl^ 


—2k + isQ —(2k — U) -\- isQ 


This is reduced to 


Qinter(w) — ^^2^ ^ ^2 j 


k^dk 


g/3/c g /3 ^ 
2k T t/ T 


e + e ^2 


+ 


g^/c _^g/32 g/3/c_^g/32 


+ 


2k + 

(e-/3/c 


(-e^^ + e-^^) 




(S43) 


2k + isQ 2k — U + isQ 
\ / 

After the analytic continuation id ^ uj -\- ie, poles that contribute to the imaginary part of the /c-integral are 


k = 


uj — U uj ujU U — uj 


2 ’ 2 ’ 2 ’ 2 

Thus the interband contribution to the optical conductivity at the finite temperature is given by 


(S44) 


(S45) 


= Im 


Winter (^) 


-lUJ 


1 

bTrcj 


(e^ 2 e ^2)(_e^ 2 g p 2 ) 

2 /IA\ —2 


CJ- t/ 


CJ-t/ 


-2 


d(u - U) 


+ ^ q oj-U fUJ\^ ^ f^\ 

+ 2(e-^^-e^^)(-) Z(-) 
+ (e-/5^+e-^t)(e 
_(e-/3-^+e-^¥)(e 


2 — e 


«)(^)p(^) 


-2 


-CJ + t/ 


-CJ + t/ 


-2 


e(-uj^U)\. (S46) 


The minus sign for the term in the last line arises because the pole k = ^ ^ locates in the lower half plane while 
other poles locate in the upper half plane. 


2. Intrahand transition 

The intraband contribution to Q(ift) is given by 

T -Sintra m], (S47) 
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where 


and 


^intra('^^) — / OdOdcj) ^ ^ G -\—|- (^5 iuJ^i^ (“l“|^a:|“l“)^H—1-(^ Q'l “1“ (“1“ |^a: | “1“) 

iujn 

_QU.2'^F{hk -^)- npihkJrq “ 


47r 1 

Yz2 L 




iH -\- hj^ hj^-^q 


J3h I ^-I3^\{^-I3h I -/3^\''^F{hk ^) npikk^q + ^) 


u 


iH H“ hj^ hj^-^q U 

+ (g/3?i + g-/?^^2 ^P(^fc + f) ~ »F(fefc+g + f) 

ill “h hj^ h/j-i-g 

+ (e/5'1 _,_ g-/3f ^ /5.F(fefc + f) -npikk+q - f) 


ill + h/c — + 1/ 


-^intra(il^) — OdOdcj) ^ ^ Cr-i^n) ( |^a:| )1^-i^n “l~ iHm) ( l^a: | ) 

iuJn 

-6^\2^F{-hk - f) - ^F(-^/c+g - f) 


47r 1 

TY L 




ill hj^ H“ h/j-i-gf 

+ (e/5'^ + e“^T)(e-/5'‘ + - f) - nF(-/ifc+g + f) 

ill hj^ “h h/j-i-gf 1/ 

+ (^g-ph ^-p^^^2 '^F-hk + f) - nF{-hk+q + f) 

ill h]^ “1“ h/j-i-g 

+ (e'®'" + e~^^)(e~^^ + -npi-hk+g - f) ^ 

ill — h]^ -|- hj^-^q “h U 


(S48) 


(S49) 


After performing an analytic continuation iH^ ^ cj+ie and taking a limit g ^ 0, we obtain the intraband contribution 
to the optical conductivity 


_intra 


(OJ) 


y’dfcfc2_^(e/5'^ + e-^¥)(e-55'^ + e-55f) 

TT\ / /"A! 

6{uj-U). (S50) 


1 

"Z7 


( ^/c - + ^F ^-^/c - - ^F ^-h/c + — 


We note that we used the equation n'^(e) = n'p{—e) in the first term, and the forth term in Eq. (S49) can be discarded 
after analytic continuation because of a factor 5{uj ^U). 

We show the temperature dependence of weights of peaks at cj = 0 and U in Fig. SI. 


D. Temperature dependence of Drude weight 

We study the behavior of the Drude weight in the limit T ^ 0. The Drude weight is given by the coefficient of 
5{uj) in Eq. (S50) as 


WDrude = ^ |-(e ^'* + e -{e^^+ e '^b^'^F | 

= [dkk^ _(S51) 

6 f 2 J (e-/ 3 fe + e/ 5 fe + 2 e-/ 5¥)2 ^ ^ 

In the noninteracting case {U = 0), the Drude weight behaves as WDrude This is obtained from a crude 

estimation by replacing the factor integrand with 1 for k < T and with 0 otherwise. On the other 
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T 


FIG. SI. Temperature dependence of weights of intraband transitions at cj = 0 (blue) and oj = U (red). 


hand, in the case of strong interactions {U oo), the Drude weight behaves as W^rude ^ e Thus the Drude 

weight is suppressed exponentially as the interaction U increases. 


E. Hall conductivity 

We study the Hall conductivity for a fixed value of kz. The Hall conductivity has a nonzero contribution from a 
combination 


{Trb+{T){+\a^\-)bl{T)b-{-\ay\+)bl). 


(S52) 


Here, we set the momentum transfer as q = 0 because we focus on the dc Hall conductivity. We note that other 
combinations of current matrices vanish. After we integrate over the direction 0 of {kx, ky) = k\\ (cos 0, sin 0) in current 
matrices [Eq. (S30)], the expectation value is given by 


Q(t) = ^^{b+{T)bl{T)b-bl) + {b-{T)bl{T)b+bl) 


k 

2'Kikz 1 
k Z 

2'Kikz 1 . 2rh+l3h , 2rh-^h\ 

k ^ ’■ 


({0\b-b+(T)bl(T)b-blbl\0)e'^^ + {0\b+b-{T)bl{T)b+blbl\0)e-'^^) 


(S53) 


With the Fourier transformation, we obtain 


Q{iujn) = dre“"^(6+(T)6L(T)6_6t) + (6-(T)6t('r)&+&t) 

2Kikz 1 /—— 

k Z \ iujn — 2/i iujn -\-2h J 

^ 2nik, 1 2iu)n 

k ^ ^-(zw„)2+4/i2 


(S54) 
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By performing analytic continuation and taking the zero frequency limit, we obtain the Hall conductivity 

aM = ^ / hdhM^) 

27r J 2k^ e^h ^-i3h 


(S55) 


where k\\ is the radial coordinate for {k^^ky) and k = ^^k\\ + k‘^. In the zero temperature limit, the Hall conductivity 
reduces to 


cr. 




oo K 

/ k^dk ^ 
0 


2fc3 


1 kz 
27r 2 k 


k\\ =oo 


k\\=0 


= ^sgn(/c,). 


If we restore the unit of the Hall conductivity is given by 

^2 


^xy 


(kz) = ^Sgn{k^), 


which remains quantized into j^h in the WML 


(S56) 


(S57) 


SIV. STABILITY OF THE MOTT GAP 


We study the stability of the Mott gap against the interaction 

iic=Y. 

k,k',q 

We consider the self-energy arising in the second order of this interaction, 

S(fc,zw)= f d^q'^V{q)V{-q)G{k + q,iuj + in)U{q,in), 

iQ 

with the density-density correlation function 

n(q, in) = d?k' Tr[G{k', ioj')G{k' - q, iuj' - zQ)]. 


(S58) 


(S59) 


(S60) 


This is explicitly written as 


E(fe,iw)= f d^qd^k' Y ^(9)^(-9) 


(zw + zf2) + (|fe + q| + ^)n{k + q) ■ cr 


X 


{iuj -h - {\k + gr| -h {ioo'y - {\k'\ -h f 


{iuj'- iny - {\k'- q\ + ^y 


X 2 


— iO) + ( \k'\ -f — ^1 + -^ ) Ti{k') • n{k' — q) 


(S61) 


If the instability for the Mott gap were present, the gap should close at /c = 0 by the consideration from the rotation 
symmetry. Therefore, we focus on the self-energy for = 0 . By summing over Matsubara frequencies and setting 
k = 0, we obtain 


'E{k = 0, iuj 



d^qd^k' 


\v{q)\^ 

2 


f, k'{k'-q) \ 


{iuj — |/c'| — \k' 


X 


iuj — \k'\ — \k' 


q\-U^ 



a 


1 

q|-[/)2-(|q| + f)2 


(S62) 


After performing an integration over k\ the terms \k' — q\ and k' • {k' — q) no longer have a dependence on the angle 
of gr, because they only depend on the relative angle between k' and q. Then the only term depending on the angle 
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of q after the k' integration is riq • cr, which vanishes upon the integration over the angle of q. Thus the self-energy 
= 0 , ioj) is diagonal with respect to the spin degrees of freedom. Furthermore, the imaginary part of uj) (after 
the analytic continuation) appears only at cj = \k'\-\- \k' — q\ \q\ ^ ^ and uj = \k'\-\-\k' — q\ — \q\ ^ 

The imaginary part of 'E{k = 0,co’) is zero for cj < Therefore, the gap of ^ in the Green’s function is stable against 
the inclusion of the interaction He- 

We note that the the perturbation theory with respect to V{q) is valid because of the absence of the infrared 
divergence. In the case of the contact quartic interaction V{q) = V, we notice that the infrared divergence does 
not appear for iuj = 0 because of the gap of ^ in the energy denominator. In the case of the repulsive Coulomb 

interaction V(q) = the infrared divergence is also absent, because the density-density correlation function 

behaves n(gr, iQ) oc for small q and II, and the integral is convergent around g = 0. 



